Gravitational force in weakly correlated particle distributions 
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Abstract 

We study the statistics of the gravitational (Newtonian) force in a particular kind of weakly corre- 
lated distribution of point-like and unitary mass particles generated by the so-called Gauss-Poisson 
point process. In particular we extend to these distributions the analysis a la Chandrasekhar in- 
troduced for purely Poisson processes. In this way we can find the asymptotic behavior of the 
probability density function of the force for large values of the field as a generalization of the Holtz- 
mark statistics. The validity of the introduced approximations is positively tested through a direct 
comparison with the analysis of the statistics of the gravitational force in numerical simulations of 
Gauss-Poisson processes. Moreover the statistics of the force felt by a particle due to only its first 
nearest neighbor is analytically a numerical studied, resulting to be the dominant contribution to 
the total force. 
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I. INTRODUCTION 

The knowledge of the statistical properties of the grav- 
itational field in a given distribution of point-particles 
with the same (unitary) mass provides an important in- 
formation about the system in many cosmological and as- 
trophysical applications in which the particles are treated 
as elementary objects. In fact such an information ac- 
quires a particular importance in the contexts of stel- 
lar dynamics and of the cosmological n— body simula- 
tions to study the formation of structures from initial 
mass density perturbations 1 . .2] . Similar studies are in- 
volved in other domains of Physics, such as the analysis 
of the statistics of the dislocation- dislocation interaction 
for what concerns the analysis of crystal defects in con- 
densed matter physics Until now a complete study 
of this problem has been accomplished only in the case 
of an uncorrelated Poisson particle distribution 0, Im- 
partial results have been found more recently in other 
two cases: 1) a fractal point distribution 0, and 2) a 
radial density of particles Q. In this paper we study 
the case of a so-called Gauss-Poisson (GP) point process 
[Til IT^ . which generates particle distributions charac- 
terized by only two-points correlations, i.e. connected 
n— points correlation functions vanish for n > 3. In this 
sense it can be seen as the first step of correlated systems 
beyond the completely uncorrelated Poisson distribution. 
For the GP point process we generalize the method used 
by Chandrasekhar [Bj] for the Poisson case introducing 
some approximations. Moreover we study the contribu- 
tion to the total force experienced by a particle due to 
the first nearest neighbor (NN), in order to evaluate the 
weight of the granular neighborhood of a fixed particle. 
Indeed, as clearly discussed by Chandrasekhar ||g, one 
the main problem of the dynamics of a self gravitating 
particle distribution is concerned with the analysis of the 



force acting on a single particle. Such a study is at the 
basis of the analysis of particle and fluid dynamics. In 
a general way, it is possible to show that there are two 
different contributions: the first is due to the system as 
a whole and the second is due to the influence of the im- 
mediate neighborhood of the particle. The former is a 
smoothly varying function of position and time while the 
latter is subject to relatively rapid fluctuations. These 
fluctuations, which are then the subject of the present 
paper, are related to the underlying statistical properties 
of particle distribution. 



II. GRAVITATIONAL FORCE PROBABILITY 
DENSITY IN A POISSON DISTRIBUTION 

Firstly, let us recall Chandrasekhar's Q results for the 
Poisson case. A Poisson distribution of point-particles 
with average density n in a volume V is obtained by oc- 
cupying randomly with a particle of unitary mass each 
volume element dV with probability ndV (with n > 
and ndV <C 1) or leaving it empty with complementary 
probability 1 — ndV with no correlation between different 
volume elements. Therefore the average number of par- 
ticles in the volume V is N — nV with fluctuations from 
realization to realization of the order of \J nV (the so- 
called Poissonian fluctuation), which become negligible 
with respect to N in the large V limit. By definition the 
connected two-point correlation function £(r) has only 
the diagonal part, that is £(r) = 8(f) /n. Any other sta- 
tistically homogeneous particle distribution is character- 
ized by a connected two-points correlation function of the 
form [111 



n 



(1) 



2 



where £(r) is the non-diagonal part due to correlations 
between the positions of different particles. 

In general, given the particle distribution, the gravita- 
tional field acting on the origin of axis is given by: 
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where the sum runs over all the system particles. Once 
the statistical ensemble of particle distributions is chosen, 
one can evaluate the probability density function (PDF) 
W(F) of the field F by taking the average of 5(F— ^ ^) 
over the ensemble. In particular for a Poisson distribu- 
tion in a volume V with average density of particles n 
this calculation can be performed in the following way 
PJ. Since in this case the positions of different particles 
are not correlated at all, the joint PDF of the positions 
of the N particles of the system is simply given by: 



P{r 1 ,r 2 , 
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Therefore the PDF of the total gravitational force acting 
on the origin is: 
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By using the Fourier representation of the Dirac delta 
function and taking the thermodynamic limit V, N — > 
+oo with N/V — n, we obtain: 



W{F) 



where 
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Note that A(k) = cxp(— nCp{k)) is the Fourier transform 
of W{F), i.e. A{k) is the so-called characteristic function 
of the stochastic field F p|. Since Cp(k) depends only 
on k = \k\, Eq. 0] says that the direction of F is 

completely random (i.e. with a PDF p(6,<fi) — 1/(4-71")) 
and decoupled from F = \F\, whose PDF (defined with 
F > 0) is instead given by 



W(F) 
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This important result is known under the name of Holtz- 
mark distribution. An explicit expression of W(F) is not 
possible to be obtained; anyway it is rather simple to 
study the asymptotic regimes. Probably the most im- 
portant feature of Eq. HJ is the behavior at large F that 
can be found to be W{F) ~ 2imF~§ for F > n~ 2/3 . 
Now we show that this limit behavior is completely de- 
termined by the position of the (NN) particle. In order 



to show this result we have to evaluate the probability 
u>(r)dr that, given a particle, its first NN particle is at a 
distance between r and r + dr. An equation satisfied by 
u)(r) can be found by considering that the probability of 
finding the first neighbor between r and r + dr is equal to 
the product of the probability that there is no particle in 
the distance interval (0, r] and the probability Amnr^dr 
of finding a g eneric particle in the interval of distances 
(r, r + dr] [l(j , that is 
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The derivation of Eq.[7|is based on the fact that there is 
no correlation between the position of different particles, 
implying that the probability of finding no particle in 
(0, r] is independent of the probability of finding a parti- 
cle in (r, r + dr] . This of course holds for a homogeneous 
Poisson distribution, but in general it is not true for cor- 
related distributions. Equation0can be simply solved to 



give: 
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By considering that F — 1/r 2 , we can find, by a sim- 
ple change of variable, the PDF of the modulus of the 
gravitational field generated by the first neighbor as: 



W nn (F) = 2TrF-inexp[ - 



In the limit F > n 2/3 Eq. El reads 
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W nn (F) ~ 2TmF- 
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which is exactly the same of the asymptotic behavior of 
the PDF W{F) of the modulus of the total force. This 
result implies that in a Poisson distribution (see Fig^l, 
in which large-scale correlations are absent but density 
fluctuations are present at all scales, the main contribu- 
tion to the force acting on a particle comes from particles 
in its neighborhood, the rest of faraway particles giving 
only a finite additional contribution because of statistical 
isotropy. 



III. THE GAUSS-POISSON POINT PROCESS 

We now discuss the statistical properties of the gravi- 
tational Newtonian field arising by an infinite weakly cor- 
related particle system, and in particular of a so-called 
Gauss-Poisson distribution of point-like field sources (of 
unitary mass). A GP particle distribution 0,0] is built 
in the following way: first of all, let us take a Poisson 
distribution of particles of average density no > 0. The 
next step is to choose randomly a fraction < q < 1 
of these Poisson points and to attach to each of them 
a new "daughter" particle in the volume element d 3 r at 



vectorial distance r from the "parent" particle with prob- 
ability p(r)d 3 r independently one each other. Therefore 
the net effect of this algorithm is of substituting a frac- 
tion q of particles of the initial Poissonian system with an 
equal number of correlated binary systems. This is the 
reason why this kind of point distribution can be very 
useful in all the physical applications characterized by 
the presence of binary systems. 

It is simple to show that the final particle density in the 
so-generated GP distribution is n = n (l + q). It is also 
possible to show that the connected two-point correlation 
function is 
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and that all the other connected n— points correlation 
function with n > 3 vanish [13|. This means that all the 
statistics of a GP stochastic distribution is reduced to the 
knowledge of one and two-points correlations. For this 
reason the GP particle distribution is the discrete anal- 
ogous of the continuous Gaussian continuous stochastic 
fields. Moreover since p{r) is a PDF, £(r) is non- negative 
and integrable over all the space, hence correlations are 
short ranged. This is the reason why the GP point- 
particle distributions can be seen as the most weakly 
correlated particle system beyond the Poisson one. To 
show the validity of Eq. ^] is a quite simple task, in fact 
it is sufficient to use the definition of average conditional 
density n p (r) of particles seen by a generic particle of the 
system at a vectorial distance r from it without counting 
the observing particle itself. It is simple to show that : 



n p (r) = n [1 + £(r)] 



(12) 



where £(r) is the non-diagonal part of On the other 
hand in the GP model the conditional average density can 
be evaluated in the following way: the number of particles 
seen in average by the chosen particle in the origin in 
the volume element d 3 r around r is nd 3 r if the chosen 
particle is neither a "parent" nor a "daughter" (i.e. with 
probability (1— q)/(l+q)) and nd 3 r+p(f)d 3 r if it is either 
a "parent" or a "daughter" (i.e. with a complementary 
probability 2q/{\ + q)). This gives directly 
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which is equivalent to Eq. ^] Note that if p(r) depends 
only on r (i.e. it is spherically symmetric) then the par- 
ticle distribution, in addition to be statistically homoge- 
neous (i.e. translational invariant), is also statistically 
isotropic (i.e. rotational invariant). 



IV. GENERALIZATION OF THE HOLTZMARK 
DISTRIBUTION TO THE GAUSS-POISSON 
CASE 
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FIG. 1: Connected two-points correlation function plus one 
measured in different types of the particles distribution. The 
correlations function are taken by averaging over 1000 real- 
izations of a sample of N = 10 2 particles (in the same fixed 
volume V — 1.) Poisson case, the Gauss-Poisson case (GP- 
box) where p(r) is a box — function, and the Gauss-Poisson 

case (GP-Exp), where p(r) 
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erated a GP distribution with fixed n > and < q < 1 
in a volume V (therefore N — nV with negligible relative 
fluctuations in different realizations of the distribution in 
the infinite volume limit). Let us also set the coordinate 
system in such a way that the origin is occupied by a 
particle of the system. We want to calculate the PDF 
W(F) of the total gravitational field F acting on the ori- 
gin of coordinates due to all the particles out of the origin 
conditioned to the fact that this point is occupied by a 
particle of the system . Therefore if the particles seen 
by the one in the origin are N — 1 and p c ifi, ^2, fjv— l) 
is the PDF of their position with respect the origin we 
can write: 
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Pc(fi,f 2 , ...,rW_i) • (13) 
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The main problem to face is due to the fact that, since 
in the GP case, two-points correlations are present, 
p c (ri, fa, r/v-i) cannot be written as a product N — 1 
one particle PDF's as in the Poisson case. This prop- 
erty would prevent from the possibility of applying the 
Markov method to this case, and an explicit evaluation 
of W would become impossible. For this reason we intro- 
duce the approximation consisting in imposing the fac- 
torization 

N 

Pc{ri,f 2 , ...,fjv_i) = JJr(fi) , 



We can now try to generalize the Holtzmark distribu- 
tion to this correlated case. Let us suppose of having gen- 



but taking into account the fact that in average the par- 
ticle in the origin sees a density of particles in the point 
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'given by Eq. ^]with £(f 
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This approximation permits to use the Markov method 
to find W(F) which, in the limit V — > +oo, can be shown 
to be given by 



IT(F) 



where 
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d 3 Fexp(i/c-F-nC GP (A:)) , (15) 
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As shown below by a direct comparison with the results 
of numerical simulations, this approximation is quite ac- 
curate at least in the asymptotic regime F — > +oo. Note 

that the function A(k) — exp (— nCGp(k)j is nothing 

else the characteristic function of the total force F act- 
ing on the particle in the origin in the GP case. As 
aforementioned, if the PDF p(r) depends only on r = \r\, 
the particle distribution is statistically isotropic. Conse- 
quently, as for the Poisson case, W(F) also depends only 
on F — \F\ and A(k) on k = \k\. That is the direction 
of F is completely random while the PDF of F is given 
by (by recalling p(r) with p{r) to put in evidence the 
dependence only on r) 



W{F) = 4ttF 2 W(F) 
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We limit the rest of the discussion to this isotropic case. 
As for the Poisson distribution, it is not possible to find 
an explicit form of W(F) (or equivalently of W(F)). 
However we can connect their large F behavior to the 
small r behavior of p(r) and to that of the Poisson case. 
In order to do this, it is important to use the general 
properties of the Taylor expansion of the characteris- 
tic function A(k) to the lowest order greater than zero. 
In particular in this isotropic case we use that if 
W(F) ~ F~ a at large F (note that a > 3 in any case as 
W(F) is a normalizable PDF) then 



A(k) = / d 3 Fexp (-ik ■ F) W{F) = 



\{F 2 )k 2 ifa>5 
ak a - 3 if 3 < a < 5 



(18) 



where a > is a constant characterizing the singularity. 
Note that a > 5 implies that (-F 2 ) is finite, and that for 



the Poisson case a — 9/2, and correspondingly A(k) ~ 

l-if(27Tfc)l. 

Therefore our strategy is to find a by connecting the 
expansion given in Eq. 1181 to the form of p(r) and in 
particular to its small r behavior. Let us suppose that 
p(r) ~ r@ at small r (in any case j3 > —3 as p{r) is a 
PDF of a three-dimensional stochastic variable). It is 
quite simple to show that at small k the integral 



l{k;0) 



behaves as follows: 
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Cyk 2 if < 1 
c 2 k 2 if p > 1 . 
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where c\ and ci are two positive constants depending on 
the specific form of p(r) . Consequently, by inserting this 
result in Eq. El we can distinguish three cases for what 
concerns the asymptotic behavior of W(F): 

• For p > the dominating part in A(k) at small k 
is exactly the same as in the Poisson distribution 
with the same average density n, i.e. 



A(k) 
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4n 
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which implies W(F) ~ F 2 (or equivalently 
W{F) = 4ttF 2 W(F) ~ F"S) at large F with the 
same amplitude of the pure Poisson case. 

• For p = we have again a substantially Poisso- 
nian behavior but the coefficient of the non zero 
order term receives a contribution from two-points 
correlations: 



A(k) ~ 1 - 87m 
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n(l + q) 
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which implies again W(F) ~ F _ s at large F but 
with a larger amplitude than in the Poisson case. 

• for p < the small k behavior of A(k) is completely 
changed, being 
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F 2 (or equivalently 



which gives W(F) 
W{F) = 4:nF 2 W(F) - F-^) 



V. COMPARISON WITH SIMULATIONS 

In order to check the validity of these theoretical re- 
sults, we have performed numerical simulations consist- 
ing in generating two kinds of GP distributions of parti- 
cles with two explicite choices of p(r) (see Fig.^l, and in 
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simulation 
theoretical 




FIG. 2: Comparison between theoretical predictions (con- 
tinuous lines) and simulations of the tail of the probability 
density function of the modulus of the gravitational force 
W(F) = 4nW(F). The average for the PDF is performed 
over 2 ■ 10 4 realizations of a the same distribution with each 
of N = 10 5 particles. Top panel: Poisson case. Middle panel: 
Gauss-Poisson case, where p(r) is a box — function. Bottom 
panel: Gauss-Poisson case, where p(r) has an exponential 

cut-off: p(r) — 



exp(- 



r) 



47TTQ 



measuring directly the PDF of W(F) in these cases: 
(1) In the first one the choice of p{r) is simply a positive 
constant up to a fixed distance ro and zero beyond this 
distance: 



p(r) 







if < r < r 
if r < ro . 



(23) 



That is the probability of attaching a "daughter" particle 
at a distance between r and r + dr from its "parent" is 
3r 2 dr/r$ if r < ro and zero for r > r$. As shown above 
this choice of p(r) should give W(F) ~ F~i at large F 
but with a larger amplitude than the pure Poisson case. 
(2) In the second case p(r) decays exponentially fast at 
large r but increases as r~ 2 at small r, i.e. 



p{r) 



1 exp(- 
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This choice of p(r) should give W(F) ~ F 2 at large 
F. The results of these simulations for the large F be- 
havior of W(F) have been compared to the theoretical 
previsions showing a good agreement (see Fig. 0). Con- 
sequently the validity of the approximation at the base of 
these calculations has been positively tested. The found 
validity of the relations between the small scale behavior 
of £(r) and the large scale behavior of W(F) for the GP 
case suggests that this approach can be extended to more 
general cases of correlated particle distributions. 



VI. NEAREST NEIGHBOR APPROXIMATION 
FOR THE GAUSS-POISSON CASE 

As for the Poisson case, we now analyze, for GP distri- 
butions in the statistical isotropic case, the importance 
of the first NN contribution to the total force felt by the 
particle in the origin. To this aim we use again only the 
information that the average conditional density n p (f) 
(which depends only on r) seen by the particle in the 
origin is given by Eq. 1121 Following exactly the same 
reasoning done for the Poisson case with n p (r) replacing 
the simple n, we can write: 
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where again £(r) 
solved to give 



u)(x)dx) inirr 2 (1 + £,(r)) , (25) 
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where we have called ujp(r) the function ui{r) found above 
for the Poisson case. By imposing p(r) ~ r@ at small r 
and using again F = 1/r 2 in order to pass from u>(r) to 
W nn (F), it is simple to see that W nn {F) has the same 
aforementioned scaling behavior at large F of W(F) for 
all the permitted values of j3 and with the same coef- 
ficient. Therefore also in the GP case we have that the 
main contribution to the force felt by a particle in the sys- 
tem is due to its first NN. This was somehow expected be- 
cause the main change introduced by passing from Pois- 
son to GP distribution concerns mainly the introduction 
of additional fluctuations of density in the neighborhood 
of any particle. 



VII. DISCUSSION 

In conclusion we have studied the gravitational force 
distribution in the so called Gauss-Poisson particle distri- 
butions, which can be considered in some sense the most 
weakly correlated case beyond the Poisson one. For this 
particle systems we have seen how to generalize the meth- 
ods developed for the Poisson case in order to find the 
PDF of the gravitational force. The main result is that, 
in the GP case, significative deviations from the Poisson 
behavior can be caused only by the small scale behavior 
of two-points correlations, which can introduce strong 
modification in the large force regime when diverging at 
small distances. This is confirmed by direct results in 
numerical simulations. Moreover, as in the Poisson case, 
we found that the main contribution to the force felt 
by a generic particle is due mainly to its neighborhood. 
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The importance of this work is twofold, (i) firstly, this 
is the first case of statistically homogeneous correlated 
particle distribution in which a systematic study of the 
gravitational force a la Chandrasekhar is done, (ii) this 
study suggest some basic ingredients to be used in fu- 
ture attempts of extending the analysis to more complex 
correlated particle distributions. 
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